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Abstract 

We present parameter-free LDA+DMFT (local density approximation + dynamical mean field 
theory) results for the many-body spectra of cubic SrVC>3 and orthorhombic CaV03. Both systems 
are found to be strongly correlated metals, but not on the verge of a metal-insulator transition. 
In spite of the considerably smaller V-O-V bond angle in CaV03 the LDA+DMFT spectra of the 
two systems for energies E < Ep are very similar, their quasiparticle parts being almost identical. 
The calculated spectrum for E > Ep shows more pronounced, albeit still small, differences. This is 
in contrast to earlier theoretical and experimental conclusions, but in good agreement with recent 
bulk-sensitive photoemission and x-ray absorption experiments. 

PACS numbers: 71.27.+a, 71.30.+h 
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I. INTRODUCTION 



Transition metal oxides are an ideal laboratory for the study of electronic correlations in 
solids. In these materials, the 3d bands are comparatively narrow with width W ~ 2 — 3 eV 
so that electronic correlations, induced by the local Coulomb interaction U m 3 — 5 eV, 
are strong. For U/W <C 1, we have a weakly correlated metal and the local density 
approximation (LDA) (see for a review) works. In the opposite limit U/W 3> 1 and for an 
integer number of 3d electrons, we have a Mott insulator with two separate Hubbard bands as 
described by Hubbard's I and III approximations j^j or the LDA+U method Transition 
metals are, however, in the "in-between" regime, U /W = 0(1), where the metallic phase is 
strongly correlated with coexisting quasiparticle peak at the Fermi level and Hubbard side 
bands. 

Dynamical mean-field theory (DMFT) jj, S, 0, Q, 0| , a modern, non-perturbative many- 



ike the one-band Hubbard model (for an 
). And, with the recent merger P.llollri| 



body approach, is able to capture the physics in the whole range of parameters from 
U/W <C 1 to U /W » 1 for model Hamiltonians, 
introduction into DMFT and its applications see 
of LDA and DMFT, we are now able to handle this kind of correlation physics within realistic 
calculations for transition metal oxides. 

A particularly simple transition metal oxide system is Cai-^Sr-rVOs with a 3d 1 electronic 
configuration and a cubic perovskite lattice structure, which is ideal for SrVOs and or- 
thorhombically distorted upon increasing the Ca-doping x. Fujimori et al. 12j initiated the 
present interest in this 3c? 1 series, reporting a pronounced lower Hubbard band in the pho- 
toemission spectra (PES) which cannot be explained by conventional LDA. Subsequently, 
thermodynamic properties were studied, and the Sommerfeld coefficient, resistivity, and 
paramagnetic susceptibility were found to be essentially independent of x (l3j . 

On the other hand, PES Q], seemed to imply dramatic differences between CaV03 and 
SrVOa, with CaVOs having (almost) no spectral weight at the Fermi energy. Rozenberg et 
al. |l5j] interpreted these results in terms of the (DMFT) Mott-Hubbard transition in the 
half-filled single-band Hubbard model. With the same theoretical ansatz, albeit different 
pa—, the optica, conductivity Q and the PES experts Q for different doping 
x were reproduced. 

The puzzling discrepancy between these spectroscopic and the thermodynamic measure- 
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ments, suggesting a Mott-Hubbard metal-insulator transition or not, were reso 
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ved recently 
. In the lat- 



by bulk- sensitive PES obtained by Maiti et al. |17[ and Sekiyama et al 
ter work it was shown that (i) the technique of preparing the sample surface (which should 
preferably be done by fracturing) is very important, and that (ii) the energy of the X-ray 
beam should be large enough to increase the photoelectron escape depth to achieve bulk- 
sensitivity Theoretically, pronounced differences between SrVC>3 surface and bulk spectra 



were reported by Liebsch |20j]. Besides bu 
strumental resolution (about 100 meV in 



k-sensitivity, the beam should provide a high in- 
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19j). With these experimental improvements, 
the PES of SrVC>3 and CaVC>3 were found to be almost identical 12, 3,0]) implying consis- 
tency of spectroscopic and thermodynamic results at last. This is also in accord with earlier 
Is x-ray absorption spectra (XAS) by Inoue et al. 2lJ which differ only slightly above the 
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Fermi energy, as well as with BIS data Q]. Parameterfree LDA+DMFT results by us 
reported in the joint experimental and theoretical Letter by Sekiyama et al. showed 
good agreement with the high-resolution bulk-sensitive PES. Independent LDA+DMFT cal- 
culations for 3d 1 vanadates, including SrVC>3 and CaVC>3 have been reported by Pavarini et 
al. j^J. Most recently Anisimov et al. employed a full-orbital computational scheme based 
on Wannier functions formalism for LDA+DMFT to calculate total densities of states for 
SrVC>3 which are in good agreement with new PES data. j^J 

In this paper we present details of our LDA+DMFT calculation j2] which could not be 
included in the joint Letter Q|. In particular, we discuss the crystal structure of SrVC>3 
and CaVC>3, explain the origin of the surprisingly small change in bandwidth of these two 
systems, and compare the calculated spectrum for energies E > Ep with XAS data [2lj . We 



start with discussing the structural differences between SrVC>3 and CaVC>3 in Section Til Al 
based on the most recent crystallographic data provided by Inoue 2f| . How these differences 



reflect in the LDA results is shown in Section III Bl before we turn to the LDA+DMFT 
scheme in Sections IIII Al and the effects of electronic correlations in Section lill Bl Section ITVl 
concludes our presentation. 
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II. SRVQ 3 VS. CAV0 3 : STRUCTURAL DIFFERENCES AND LDA RESULTS 



Structural differences 

Both SrV03 and CaV03 are perovskites, with SrV03 in the ideal cubic structure 
(Fig. [IJ and CaV03 in the orthorhombicaly distorted (GdFeOs) Pbnm struc- 
(Fig. |2J). The ideal cubic structure of SrV0 3 implies that the V0 6 octahedra, 



Pm3m 
ture 
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the main structural elements of the crystal, are not distorted: The V-Obasai and V-O ape x 
distances are the same, the O apex -V-Ob a sai and Obasai-V-Obasai angles are 90°, and the V-O-V 
bond angles are exactly 180°. The substitution of the large Sr 2+ ions by smaller Ca 2+ ions 
leads to a tilting, rotation, and distortion of the VOg octahedra. Nonetheless, in CaVC>3, 
both the V-Obasai and the V-O ap ex distances are practically the same. Therefore, the most 
important feature of the distortion is the change of the V-Ob asa i-V and V-O ape x-V angles, 
which determines the strength of the effective V3d-V3d hybridization and, hence, the band- 
width. Due to the rotation and tilting of the octahedra both angles (one of them is marked 
as Z123 in Figs. Q and 12} have the same value 9 = 180° for SrV03 and Obasai ~ 161°, 
P iane « 163°for CaV0 3 . 

B. LDA results 

The valence states of SrV0 3 and CaV0 3 consist of completely occupied oxygen 2p bands 
and partially occupied V 3c? bands. Since the V ion has a formal oxidation of 4+, there is 
one electron per V ion occupying the ^-states (configuration d l ). In Figs. El and 0] the LDA 
density of states (DOS) are presented for both compounds, as calculated by the linearized 
muffin-tin orbitals method (LMTO) 2^| with an orthogonal basis set of Sr(5s, bp, 4d, 4/), 
Ca(4s, Ap, 3d), V(4s, Ap, 3d), and 0(3s,2p, 3d) orbitals |3j|. These results are in agreement 
with previous LDA calculations in the basis of augmented plane waves (APW) by Tagekahara 
(3^. For a combination of the LDA band structure with many-body techniques P| however, 
the atomic-like LMTO wave functions employed here are particularly convenient. As is 
apparent from Figs. OjandlH the 0-2p states lie between -7.5eV and -2.0eV in SrV03 and 
between -6.5eV and -1.5eV in CaV03. The 3c? states of V are located between -l.leV and 
6.5eV in SrV03 and between -l.OeV and 6.5eV in CaV03. These results are in agreement 
with previous LDA calculations in the basis of augmented plane waves (APW) by Tagekahara 



all, 



3 If . but for our purposes LMTO is more appropriate. 

In both compounds, the V ions are in a octahedral coordination of oxygen ions. The 
octahedral crystal field splits the V-3d states into three t2 g orbitals and two e g orbitals. In 
the cubic symmetry of SrVC>3, hybridization between t2 9 and e g states is forbidden, and 
the orbitals within both subbands are degenerate. In contrast, the distorted orthorhombic 
structure of CaVC>3 allows them to mix and the degeneracy is lifted. In the lower part of 
Figs. El and we present these t2 g and e g subbands of the V-3d band. For both Sr and 
Ca compounds, the reader can see an admixture of these Vanadium t2 9 and e g states to 
the oxygen 2p states in the energy region [-8eV,-2eV]. This is due to hybridization and 
amounts to 12% and 15% of the total t2 9 weight for SrVC>3 and CaVC>3, respectively. The 
main t 2g (e 9 ) weight lies in the energy region [-1.1 eV,1.5 eV] ([0.0eV,6.5eV]) in SrV0 3 and 
[-1.0 eV,1.5 eV] ([0.2 eV,6.5 eV]) in CaVC>3. Since the energy difference between the centers 
of gravities of the t 2g and e g subbands is comparable with the bandwidth, t 2g and e g states 
can be considered as sufficiently separated in both compounds. 

Let us now concentrate on the t2 9 orbitals crossing the Fermi energy and compare the 
LDA t2 3 DOS of SrVC>3 and CaV03. Most importantly, the one-electron t2 g LDA bandwidth 
of CaV03, defined as the energy interval where the t2 9 DOS in Figs. El and 0] is non-zero, is 
found to be only 4% smaller than that of SrVO,3 (Hcavo, = 2.5 eV, WsrV0 3 = 2.6 eV); in 
general agreement with other LDA calculations [23 1. This is in contrast to the expectation 
that the strong lattice distortion with a decrease of the V-O-V bond angle from 180° to 162° 
affects the t2 9 bandwidth much more strongly. Such a larger effect indeed occurs in the e g 
bands for which the bandwidth is reduced by 10%. To physically understand the smallness 
of the narrowing of the t2 9 bands we have investigated a larger Hamiltonian consisting of e g , 
t2 g , and oxygen p orbitals. We calculated the overall direct d-d and the indirect d^p-d hopping 
for this larger Hamiltonian. The predominant contribution to the e g -e g hopping is through 
a d—p—d hybridization, which is considerably decreasing with the lattice distortion. This is 
also the case for the t2 g orbitals. However, for the t2 g orbitals, the direct d—d hybridization 
is also important. This hybridization increases with the distortion since the t2 g orbital lobes 
point more directly towards each other in the distorted structure. Hopping parameters for 
individual orbitals change rather substantially from SrV03 to CaV03. But, altogether, a 
decreasing d—p—d hybridization and an increasing d—d hybridization results in a very small 
change of the t 2g bandwidth. 
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III. CORRELATION EFFECTS IN SRVO3 AND CAVO3 



A. Dynamical Mean-Field Theory 

As is well-known, LDA does not treat the effects of strong local Coulomb correlations ad- 
equately. To overcome this drawback, we use LDA+DMFT as a non-perturbative approach 
to study strongly correlated systems It combines the strength of the LDA in 

describing weakly correlated electrons in the s- and p-orbitals, with the DMFT treatment 
of the dynamics due to local Coulomb interactions. In the present paper we will discuss the 
relevant parts of the LDA+DMFT approach only briefly, refering the reader to Ref. lol 
for details. 

The LDA+DMFT Hamiltonian can be written as 

(1) 

mi i mj^m' aa' 

Here, the index % enumerates the V sites, m the individual X> 2g orbitals and a the spin, -^leia 
is the one-particle Hamiltonian generated from the LDA band structure with an averaged 
Coulomb interaction subtracted to account for double counting U and U' denote the 
local intra-orbital and inter-orbital Coulomb repulsions, and J is the exchange interaction. 

We calculated these interaction strengths by means of the constrained LDA method |33| 
for SrV0 3 , allowing the e g states to participate in screening |33J. The resul ting value of 
the averaged Coulomb interaction is U — 3.55 eV (U = U' for t 2g orbitals lOj, |34|) and 
J = 1.0 eV. The intra-orbital Coulomb repulsion U is then fixed by rotational invariance to 
U — U' + 2 J — 5.55 eV. We did not calculate U for CaV0 3 because the standard procedure 
to calculate the Coulomb interaction parameter between two t 2 g electrons, screened by e g 
states, is not applicable for the distorted crystal structure where the e g and t 2g orbitals are 
not separated by symmetry. On the other hand, it is well-known that the change of the 
local Coulomb interaction is typically much smaller than the change in the DOS, which was 
found to depend only very weakly on the bond angle. That means that U for CaV03 should 
be nearly the same as for SrV03. Therefore we used U = 3.55 eV and J = 1.0 eV for 
both SrVQq and CaV03. This is also in agreement with previous calculations of vanadium 
systems and experiments jlfi| . 

DMFT maps the lattice problem Eq. onto a self-consistent auxiliary impurity prob- 
lem, which is here solved numerically by the quantum Monte-Carlo (QMC) technique j^. 
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Combined with the maximum entropy method 



spectral functions (37, 



38, 



39 



this technique allows us to calculate 



A computationally important simplification is due to the fact that, in ideal cubic per- 
ovskites, the (degenerate) t 2g states do not mix with the (degenerate) e g states. In this 
particular case, the self-energy matrix £(z) is diagonal with respect to the orbital indices. 
Under this condition, the Green functions G m (z) of the lattice problem can be expressed in 
the DMFT self-consistency equation by the Hilbert transform of the non-interacting DOS 



«-w - J *d3f^ ■ p> 

This procedure avoids the rather cumbersome and problematic /c-integration over the Bril- 
louin zone by the analytical tetrahedron method |4l|. We obtain iVm (e) for SrVOs and 
CaV0 3 from the t 2s -P r ojected DOS of Figs. El and |U by truncating the contribution in the 
oxygen range below -1.5 eV and multiplying by a renormalization factor (to have maximally 
one electron per site and orbital). This procedure is employed as it resembles the three 
effective bands of t2 g character which cross the Fermi energy and are, hence, responsible for 
the low-energy physics. 



B. LDA+DMFT(QMC) results and discussion 

The calculated LDA+DMFT(QMC) spectra for SrV0 3 (right panel) and CaV0 3 (left 
panel) are presented in Fig. E] for different temperatures. Due to genuine correlation effects, 
a lower Hubbard band at about —1.5 eV and an upper Hubbard band at about 2.5 eV is 
formed, with a well pronounced quasiparticle peak at the Fermi energy. The 4% difference 
in the LDA bandwidth between SrV0 3 and CaV0 3 is only reflected in some additional 
transfer of spectral weight from the quasiparticle peak to the Hubbard bands, and minor 
differences in the positions of the Hubbard bands. Clearly, the two systems are not on the 
verge of a Mott-Hubbard metal- insulator transition. Both, SrV0 3 and CaV0 3 , are strongly 
correlated metals, though SrV0 3 is slightly less correlated than CaV0 3 with somewhat more 
quasiparticle weight, in accord with their different LDA bandwidth. As one can see from 
Fig. El the effect of temperature on the spectrum is small for T < 700 K. 

In the left panel of Fig. El we compare our LDA+DMFT spectra (300K), which were mul- 
tiplied with the Fermi function at the experimental temperature (20 K) and Gauss broadened 
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with the experimental resolution of 0.1 eV j3] with the experimental PES data obtained by 
subtracting the experimentally estimated surface and oxygen contributions. Due to the high 
photon energy (several hundred eV) ^j], the PES transition matrix elements will have no 
strong energy dependence in the energy interval of a few eV studied here, which justifyies 
that we do not take such matrix elements into account. 

The quasiparticle peaks in theory and experiment are seen to be in very good agreement. 
In particular, their heights and widths are almost identical for both SrVC>3 and CaVC>3. 
The difference in the positions of the lower Hubbard bands may be partly due to (i) the 
subtraction of the (estimated) oxygen contribution which might also remove some 3d spectral 
weight below —2 eV, and (ii) uncertainties in the ab-initio calculation of U. 

Since inverse photoemission spectroscopy (IPES) data for SrVC>3 and CaVC>3 are not yet 
available we compare j^J the calculated spectrum for energies E > Ep with XAS data by 



Inoue et al. |2l| (right panel of Fig. |B }. We consider core- hole life time effects by Lorentz 
broadening the spectrum with 0.2 eV 42|, multiply with the inverse Fermi function (80K), 
and apply Gauss broadening given by the experimental resolution of 0.36 eV 2f3]. Again, the 
overall agreement of the weights and positions of the quasiparticle and upper t 2g Hubbard 
band is good, including the tendencies when going from SrVC>3 to CaVO^. For CaVC>3 
(Cao.9Sro.1VO3 in the experiment), the LDA+DMFT quasiparticle weight is somewhat lower 
than in experiment. 

In contrast to one-band Hubbard model calculations, our material specific results repro- 
duce the strong asymmetry around the Fermi energy w.r.t. weights and bandwidths. Our 
results also give a different interpretation of the XAS than in |21| where the maximum at 
about 2.5 eV was attributed solely to the e g band. We find that also the t 2g upper Hubbard 
band contributes in this energy range. 

The slight differences in the quasiparticle peaks (see Fig. EJ) lead to different effective 
masses, namely m* /mo = 2.1 for SrV03 and m*/mo = 2.4 for CaV03. These theoretical 
values agree with m*/mo = 2 — 3 for SrVO .3 a nd CaV03 as obtained from de Haas- van 
Alphen experiments and thermodynamics [13I l^. We note that the effective mass of CaV03 
obtained from optical experiments is somewhat larger, i.e., m*/m = 3.9 
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IV. CONCLUSION 



In summary, we investigated the spectral properties of the correlated 3d 1 systems SrVOs 
and CaVC>3 within the LDA+DMFT(QMC) approach. Constrained LDA was used to de- 
termine the average Coulomb interaction as U = 3.55 eV and the exchange coupling as 
J— 1.0 eV. With this input we calculated the spectra of the two systems in a parameter-free 
way. Both systems are found to be strongly correlated metals, with a transfer of the ma- 
jority of the spectral weight from the quasiparticle peak to the incoherent upper and lower 
Hubbard bands. The calculated DMFT spectra of SrV03 and CaV03, and in particular the 
quasiparticle parts, are found to be very similar, differences being slightly more pronounced 
at energies above the Fermi energy. Our calculated spectra agree very well with recent 
bulk-sensitive high-resolution PES |l8( and XAS 21 1 data, i.e., with the experimental spec- 



trum below and above the Fermi energy. Both compounds are similarly strongly correlated 
metals; CaVOa is not on the verge of a Mott-Hubbard transition. 

Our results are in contrast to previous theories and the expectation that the strong 
lattice distortion leads to a strong narrowing of the CaVOs bandwidth and, hence, to much 
stronger correlation effects in CaVOs. While the e g bands indeed narrow considerably, the 
competition between decreasing d—p—d and increasing d—d hybridization leads to a rather 
insignificant narrowing of the t2 9 bands at the Fermi energy. This explains why the k- 
integrated spectral functions of CaV03 and SrV03 are so similar. With our theoretical 
results confirming the new PES and XAS experiments, we conclude that the insulating-like 
behavior observed in earlier PES and BIS experiments on CaVOs must be attributed to 
surface effects 
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FIG. 1: SrV0 3 , an ideal cubic perovskite, with V-O-V angle #=Z123=180°; V: large ions, Obasal : 
small bright ions, O apcx : small dark ions. 




FIG. 2: CaVOs, an orthorhombicaly distorted perovskite, with V-O-V angles Obasai ~ 161°, p i ane m 
163°; V: large ions, Ob aS ai : small bright ions, O apex : small dark ions. The local c-axis, later used 
for the definition of the orbitals, is directed along the tilted V-O apcx direction. 
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FIG. 3: DOS of SrV0 3 as calculated by LDA(LMTO). Upper panel: V-3d (full line) and 0-2p 
(dashed line) DOS; lower panel: partial DOS of V-3d fo g ) (full line) and V-3d(e 9 ) (dashed line) 
orbitals. The Fermi level corresponds to OeV. 




FIG. 4: DOS of CaV0 3 as calculated by LDA(LMTO). Upper panel: V-3d (full line) and 0-2p 
(dashed line) DOS; lower panel: partial DOS of V-3d fo g ) (full line) and V-3d(e g ) (dashed line) 
orbitals. The Fermi level corresponds to OeV. 
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Abstract 

We present parameter-free LDA+DMFT (local density approximation + dynamical mean field 
theory) results for the many-body spectra of cubic SrV03 and orthorhombic CaVOs. Both systems 
are found to be strongly correlated metals, but not on the verge of a metal-insulator transition. 
In spite of the considerably smaller V-O-V bond angle in CaVOs the LDA+DMFT spectra of the 
two systems for energies E < Ep are very similar, their quasiparticle parts being almost identical. 
The calculated spectrum for E > Ep shows more pronounced, albeit still small, differences. This is 
in contrast to earlier theoretical and experimental conclusions, but in good agreement with recent 
bulk-sensitive photoemission and x-ray absorption experiments. 
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I. INTRODUCTION 



Transition metal oxides are an ideal laboratory for the study of electronic correlations in 
solids. In these materials, the 3d bands are comparatively narrow with width W ~ 2 — 3 eV 
so that electronic correlations, induced by the local Coulomb interaction U ~ 3 — 5 eV, 
are strong. For U jW <C 1, we have a weakly correlated metal and the local density 
approximation (LDA) (see 1 J for a review) works. In the opposite limit U jW S» 1 and for an 
integer number of 3d electrons, we have a Mott insulator with two separate Hubbard bands as 
described by Hubbard's I and III approximations 2j or the LDA+U method 3j. Transition 
metals are, however, in the "in-between" regime, U /W = 0(1) } where the metallic phase is 
strongly correlated with coexisting quasiparticle peak at the Fermi level and Hubbard side 
bands. 

Dynamical mean-held theory (DMFT) DDDDD , a modern, non-perturbative many- 
body approach, is able to capture the physics in the whole range of parameters from 
U /W -v.. 1 to/ /W 1 for model Hamiltonians^like the one-band Hubbard model ^or ;m 
introduction into DMFT and its applications see ||)- And, with the recent merger tllUllUl 
of LDA and DMFT, we are now able to handle this kind of correlation physics within realistic 
calculations for transition metal oxides. 

A particularly simple transition metal oxide system is Cai^S^VC^ with a 3d 1 electronic 
configuration and a cubic perovskite lattice structure, which is ideal for SrVOs and or- 
thorhombically distorted upon increasing the Ca-doping x. Fujimori et al. initiated the 
present interest in this 3d 1 series, reporting a pronounced lower Hubbard band in the pho- 
toemission spectra (PES) which cannot be explained by conventional LDA. Subsequently, 
thermodynamic properties were studied, and the Sommerfeld coefficient .resistivity, and 
paramagnetic susceptibility were found to be essentially independent of x 

On the other hand, PES seemed to imply dramatic differences between CaVOs and 
SrVCK, with CaVOs having (almost) no spectral weight at the Fermi energy. Rozenberg et 
al. jl^J interpreted these results in terms of the (DMFT) Mott-Hubbard transition in the 
half-filled single-band Hubbard mod el. With the same theoretical ansatz. albeit different 
I op*. co ndu c tmty an d t he PES e X pe„ m en ts *» d ifa en t d op ingS 
x were reproduced. 

The puzzling discrepancy between these spectroscopic and the thermodynamic measure- 
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ments, suggesting a Mott-Hubbard metal-insulator transition or not, wereresolved recently 
by bulk-sensitive PES obtained by Maiti et al. 17j and Sekiyama et al. . In the lat- 

ter work it was shown that (i) the technique of preparing the sample surface (which should 
preferably be done by fracturing) is very important, and that (ii) the energy of the X-ray 
beam should be large enough to increase the photoelectron escape depth to achieve bulk- 
sensitivity. Theoretically, pronounced differences between SrVC>3 surface and bulk spectra 
were reported by Liebsch g£{j. Besides buhVsensitivity, the beam should provide a high in- 
strumental resolution (about 100 meV in . With these ejmermiental improvements, 
the PES of SrVC>3 and CaVC>3 were found to be almost identical |^|^,|^|, implying consis- 
tency of spectroscopic and thermodynamic results atlast. This is also in accord with earlier 
Is x-ray absorption spectra (XAS) by Inoue et al. ^] which differ only slightly above the 
Fermi energy, as well as with BIS data Parameterfree LDA+DMFT resultsby us ^^1, 
reported in the joint experimental and theoretical Letter by Sekiyama et al. showed 
good agreement with the high-resolution bulk-sensitive PES. Independent LDA+DMFT cal- 
culations for 3d 1 vanadates, including SrVC>3 and CaVC>3 have been reported by Pavarini et 
al. Most recently Anisimov et al. employed a full-orbital computational scheme based 
on Wannier functions formalism for LDA+DMFT to calculate total densities of states for 
SrVC>3 which are in good agreement with new PES data. 

In this paper we present devils of our LDA+DMFT calculation which could not be 
included in the joint Letter In particular, we discuss the crystal structure of SrVC>3 

and CaVC>3, explain the origin of the surprisingly small change in bandwidth of these two 
systems, and compare the calculated spectrum for energies E > Ep with XAS data We 
start with discussing the structural differences between SrVC>3 and CaVC>3 in Section 
based on the most recent crystallographic data provided by Inoue 2fiJ. How these differences 
reflect in the LDA results is shown in Section before we turn to the LDA+DMFT 

scheme in Sections^^^and the effects of electronic correlations in Section^^^ Section^^ 
concludes our presentation. 
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II. SRV0 3 VS. CAV0 3 : STRUCTURAL DIFFERENCES AND LDA RESULTS 



A. Structural differences 



BotliSrVOs and CaVOs are perovskites, with SrVOs in the ideal cubic structure 
Pm3mj2jl (Fig. Q| and CaVOs in the orthorhombicaly distorted (GdFeOs) Pbnm struc- 
ture (Fig. 01. The ideal cubic structure of SrVOs implies that the VC>6 octahedra, 
the main structural elements of the crystal, are not distorted: The V-Obasai and V-O ape x 
distances are the same, the O ap ex-V-Obasai and Obasai-V-Obasai angles are 90°, and the V-O-V 
bond angles are exactly 180°. The substitution of the large Sr 2+ ions by smaller Ca 2+ ions 
leads to a tilting, rotation, and distortion of the VC>6 octahedra. Nonetheless, in CaVOs, 
both the V-Obasai and the V-O a p ex distances are practically the same. Therefore, the most 
important feature of the distortion is the change of the V-Obasai-V and V-O a p ex -V angles, 
which determines the strength of the effective V3o?-V3c/ hybridization and, hence, the band- 
width. Due to the rotation and tilting of the octahedra both angles (one of them is marked 
as Z123 in Figs. Q and have the same value 6 = 180° for SrVOs and Obasai ~ 161°, 
v iane ~ 163°for CaV0 3 . 

B. LDA results 

The valence states of SrVOs and CaVOs consist of completely occupied oxygen 2p bands 
and partially occupied V 3d bands. Since the V ion has a formal oxidation of 4+, there is 
one electron per V ion occupying the (/-states (configuration d 1 ). In Figs. 0and0the LDA 
density of states (DOS) are presentedfor both compounds, as calculated by the linearized 
muffin-tin orbitals method (LMTO) with an orthogonal basis set of Sr(5s, 5p, 4d, 4/), 
Ca(4s, 4p, 3cf), V(4s, 4p, 3d), and 0(3s,2p, 3o?) orbitals These results are in agreement 

wrih previous LDA calculations in the basis of augmented plane waves (APW) by Tagekahara 
For a combination of the LDA band structure with many-body techniques 9] however, 
the atomic-like LMTO wave functions employed here are particularly convenient. As is 
apparent from Figs. Q and Q the 0-2p states lie between -7.5eV and -2.0eV in SrVOs and 
between -6.5eV and -1.5eV in CaVOs. The 3d states of V are located between -l.leV and 
6.5eV in SrVOs and between -l.OeV and 6.5eV in CaVOs. These results are in agreement 
with previous LDA calculations in the basis of augmented plane waves (APW) by Tagekahara 
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, but for our purposes LMTO is more appropriate. 
In both compounds, the V ions are in a octahedral coordination of oxygen ions. The 
octahedral crystal held splits the V-3o? states into three t 2g orbitals and two e g orbitals. In 
the cubic symmetry of SrVC>3, hybridization between t 2g and e g states is forbidden, and 
the orbitals within both subbands are degenerate. In contrast, the distorted orthorhombic 
structure of CaVC>3 allows them to mix and the degeneracy is lifted. In the lower part of 
Figs. Q and Q we present these t 2g and e g subbands of the V-3o? band. For both Sr and 
Ca compounds, the reader can see an admixture of these Vanadium t 2g and e g states to 
the oxygen 2p states in the energy region [-8eV,-2eV]. This is due to hybridization and 
amounts to 12% and 15% of the total t 2g weight for SrVC>3 and CaVC>3, respectively. The 
main t 2g (e g ) weight lies in the energy region [-1.1 eV,1.5 eV] ([0.0 eV,6.5 eV]) in SrVOs and 
[-1.0 eV,1.5 eV] ([0.2 eV,6.5 eV]) in CaVC>3. Since the energy difference between the centers 
of gravities of the t 2g and e g subbands is comparable with the bandwidth, t 2g and e g states 
can be considered as sufficiently separated in both compounds. 

Let us now concentrate on the t 2g orbitals crossing the Fermi energy and compare the 
LDA t 2g DOS of SrVC>3 and CaVC>3. Most importantly, the one-electron t 2g LDA bandwidth 
of CaVC>3, defined as the energy interval where the t 2g DOS in Figs.0and0is non-zero, is 
found to be only 4% smaller than that of SrVO^nTcaV0 3 = 2.5 eV, Ws r vo 3 — 2.6 eV); in 
general agreement with other LDA calculations 23j. This is in contrast to the expectation 
that the strong lattice distortion with a decrease of the V-O-V bond angle from 180° to 162° 
affects the t 2g bandwidth much more strongly. Such a larger effect indeed occurs in the e g 
bands for which the bandwidth is reduced by 10%. To physically understand the smallness 
of the narrowing of the t 2g bands we have investigated a larger Hamiltonian consisting of e 3 , 
t2p, and oxygen p orbitals. We calculated the overall direct d-d and the indirect d- j p-d hopping 
for this larger Hamiltonian. The predominant contribution to the e g -e g hopping is through 
a d—p—d hybridization, which is considerably decreasing with the lattice distortion. This is 
also the case for the t 2g orbitals. However, for the t 2g orbitals, the direct d—d hybridization 
is also important. This hybridization increases with the distortion since the t 2g orbital lobes 
point more directly towards each other in the distorted structure. Hopping parameters for 
individual orbitals change rather substantially from SrVOs to CaVOs. But, altogether, a 
decreasing d—p—d hybridization and an increasing d—d hybridization results in a very small 
change of the t 2g bandwidth. 
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III. CORRELATION EFFECTS IN SRVO3 AND CAVO3 



A. Dynamical Mean-Field Theory 

As is well-known, LDA does not treat the effects of strong local Coulomb correlations ad- 
equately. To overcome this drawback, we use LDA+DMFT as a non-perturbative approach 
to study strongly correlated systems 

BE30 

It combines the strength of the LDA in 
describing weakly correlated electrons in the s- and p-orbitals, with the DMFT treatment 
of the dynamics due to local Coulomb interactions. In the present paper we will discuss ^ie 
relevant parts of the LDA+DMFT approach only briefly, refering the reader to Ref. 
for details. 

The LDA+DMFT Hamiltonian can be written as 

(1) 

m i i m^m r crcr' 

Here, the index i enumerates the V sites, m the individual t 2g orbitals and a the spin. H® DA 
is the one-particle Hamiltonian generated from the LDA band structure with an averaged 
Coulomb interaction subtracted to account for double counting flj; U and U' denote the 
local intra-orbital and inter-orbital Coulomb repulsions, and J is the exchange interaction. 

We calculated these interaction strengths by means of the constrained LDA method 
for SrVOs, allowing the e g states to participate in screening The resulting value of 

the averaged Coulomb interaction is / 3.55 eV (U = V for t 2g orbitals lUl|3j|) and 
J = 1.0 eV. The intra-orbital Coulomb repulsion U is then fixed by rotational invariance to 
U I ' • '1.1 5.55 eV. We did not calculate U for CaVOs because the standard procedure 
to calculate the Coulomb interaction parameter between two t 2g electrons, screened by e g 
states, is not applicable for the distorted crystal structure where the e g and t 2g orbitals are 
not separated by symmetry. On the other hand, it is well-known that the change of the 
local Coulomb interaction is typically much smaller than the change in the DOS, which was 
found to depend only very weakly on the bond angle. That means that U for CaVOs should 
be nearly the same as for SrVOs. Therefore we used U = 3.55 eV and J = 1.0 eV for 
both SrVOa and CaVOs. This .\s also in agreement with previous calculations of vanadium 
systems and experiments 

DMFT maps the lattice problem Eq. ^ onto a self-consistent auxiliary impurity prob- 
lem, which is here solved numerically by the quantum Monte-Carlo (QMC) technique ^^|. 
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Combined with themaxinromentropy method this technique allows us to calculate 

spectral functions QQHQ- 

A computationally important simplification is due to the fact that, in ideal cubic per- 
ovskites, the (degenerate) t 2g states do not mix with the (degenerate) e g states. In this 
particular case, the self-energy matrix £(2) is diagonal with respect to the orbital indices. 
Under this condition, the Green functions G m (z) of the lattice problem can be expressed in 
the DMFT self-consistency equation by the Hilbert transform of the non-interacting DOS 



nL 0) (uj) 



G m (z) = I de- 



z - H m (z) - e 

This procedure avoids the rather cumbersome and problematic ^-integration over the Bril- 
louin zone by the analytical tetrahedron method ^Jj. We obtain Nm\t) for SrVOs and 
CaVOs from the t 23 -projected DOS of Figs.0and0by truncating the contribution in the 
oxygen range below -1.5 eV and multiplying by a renormalization factor (to have maximally 
one electron per site and orbital). This procedure is employed as it resembles the three 
effective bands of t 2g character which cross the Fermi energy and are, hence, responsible for 
the low-energy physics. 



B. LDA+DMFT(QMC) results and discussion 

The calculated LDA+DMFT(QMC) spectra for SrV0 3 (right panel) and CaV0 3 (left 
panel) are presented in Fig.0for different temperatures. Due to genuine correlation effects, 
a lower Hubbard band at about —1.5 eV and an upper Hubbard band at about 2.5 eV is 
formed, with a well pronounced quasiparticle peak at the Fermi energy. The 4% difference 
in the LDA bandwidth between SrVOs and CaVOs is only reflected in some additional 
transfer of spectral weight from the quasiparticle peak to the Hubbard bands, and minor 
differences in the positions of the Hubbard bands. Clearly, the two systems are not on the 
verge of a Mott-Hubbard metal-insulator transition. Both, SrVOs and CaVOs, are strongly 
correlated metals, though SrVOs is slightly less correlated than CaVOs with somewhat more 
quasiparticle weight, in accord with their different LDA bandwidth. As one can see from 
Fig-0 the effect of temperature on the spectrum is small for T < 700 K. 

In the left panel of Fig.0 we compare our LDA+DMFT spectra (300K), which were mul- 
tiplied with the Fermi function at the experimental temperature (20 K) and Gauss broadened 
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with the experimental resolution of 0.1 eV 1 8 J with the experimental PES data obtained by 
subtracting the experimentally estimated surface and oxygen contributions. Due to the high 
photon energy (several hundred eV) the PES transition matrix elements will have no 
strong energy dependence in the energy interval of a few eV studied here, which justifyies 
that we do not take such matrix elements into account. 

The quasiparticle peaks in theory and experiment are seen to be in very good agreement. 
In particular, their heights and widths are almost identical for both SrVOs and CaVC>3. 
The difference in the positions of the lower Hubbard bands may be partly due to (i) the 
subtraction of the (estimated) oxygen contribution which might also remove some 3d spectral 
weight below —2 eV, and (ii) uncertainties in the ab-initio calculation of U . 

Since inverse photoemission spectroscopy (IPES) data for SrVC>3 and CaVC>3 are not yet 
available we compare the calculated spectrum for energies E > Ep with XAS data by 
Inoue et al. 21 J (right panel of Fig. 0. We consider core-hole life time effects by Lorentz 
broadening the spectrum with 0.2 eV 42 J , multiply with the inverse Fermi function (80K), 
and apply Gauss broadening given by the experimental resolution of 0.36 eV 2fiJ. Again, the 
overall agreement of the weights and positions of the quasiparticle and upper t 2g Hubbard 
band is good, including the tendencies when going from SrVOs to CaVOs. For CaVOs 
(Cao.9Sro.1VO3 in the experiment), the LDA+DMFT quasiparticle weight is somewhat lower 
than in experiment. 

In contrast to one-band Hubbard model calculations, our material specific results repro- 
duce the strong asymmetry around the Fermi energy w.r.t. wekdits and bandwidths. Our 
results also give a different interpretation of the XAS than in |2j| where the maximum at 
about 2.5 eV was attributed solely to the e g band. We find that also the t 2g upper Hubbard 
band contributes in this energy range. 

The slight differences in the quasiparticle peaks (see Fig. Q lead to different effective 
masses, namely m*/m = 2.1 for SrVOs and m*/m = 2.4 for CaVOs. These theoretical 
values agree with m*/m = 2 — 3 for SrVQa and CaVOs as obtained from de Haas-van 
Alphen experiments and thermodynamics [^1^3 • We note that the effecth^mass of CaVOs 
obtained from optical experiments is somewhat larger, i.e., m*/m = 3.9 
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IV. CONCLUSION 



In summary, we investigated the spectral properties of the correlated 3d 1 systems SrVOs 
and CaVC>3 within the LDA+DMFT(QMC) approach. Constrained LDA was used to de- 
termine the average Coulomb interaction as U = 3.55 eV and the exchange coupling as 
J = 1.0 eV. With this input we calculated the spectra of the two systems in a parameter-free 
way. Both systems are found to be strongly correlated metals, with a transfer of the ma- 
jority of the spectral weight from the quasiparticle peak to the incoherent upper and lower 
Hubbard bands. The calculated DMFT spectra of SrVOs and CaVOs, and in particular the 
quasiparticle parts, are found to be very similar, differences being slightly more pronounced 
at energies above the Fermi energv. Our calculated spectra agree very well with recent 
bulk-sensitive high-resolution PES JFJ and XAS ^] data, i.e., with the experimental spec- 
trum below and above the Fermi energy. Both compounds are similarly strongly correlated 
metals; CaVOs is not on the verge of a Mott-Hubbard transition. 

Our results are in contrast to previous theories and the expectation that the strong 
lattice distortion leads to a strong narrowing of the CaVOs bandwidth and, hence, to much 
stronger correlation effects in CaVOs. While the e g bands indeed narrow considerably, the 
competition between decreasing d—p—d and increasing d—d hybridization leads to a rather 
insignificant narrowing of the t 2g bands at the Fermi energy. This explains why the k- 
integrated spectral functions of CaVOs and SrVOs are so similar. With our theoretical 
results confirming the new PES and XAS experiments, we conclude that the insulating-like 
behavior observed in earlier PES and BIS experiments on CaVOs must be attributed to 
surface effects 2(1 1. 
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FIG. 1: SrV0"3, an ideal cubic perovskite, with V-O-V angle #=Z123=180°; V: large ions, Obasal : 
small bright ions, O ap e X : small dark ions. 




FIG. 2: CaVO"3, an orthorhombicaly distorted perovskite, with V-O-V angles 9basal ~ 161°, p i ane fa 
163°; V: large ions, Obasal : small bright ions, O ap e X : small dark ions. The local c-axis, later used 
for the definition of the orbitals, is directed along the tilted V-O ape x direction. 



13 



> 

CD 

s 

o 

CO 
U 

3 



O 
Q 




Energy (eV) 

FIG. 3: DOS of SrV0 3 as calculated by LDA(LMTO). Upper panel: V-3d (full line) and 0-2p 
(dashed line) DOS; lower panel: partial DOS of Y-3d {p2g} (full line) and V-3d(e 3 ) (dashed line) 
orbitals. The Fermi level corresponds to eV. 




FIG. 4: DOS of CaV0 3 as calculated by LDA(LMTO). Upper panel: V-3d (full line) and 0-2p 
(dashed line) DOS; lower panel: partial DOS of Y-3d (i2§) (full line) and V-3d(e 3 ) (dashed line) 
orbitals. The Fermi level corresponds to eV. 
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FIG. 5: LDA+DMFT(QMC) spectrum of SrV0 3 (right panel) and CaV0 3 (left panel) calculated 
at T^300 K (solid lines), T^700 K (dashed lines), and T^llOO K) (doted lines). 
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FIG. 6: Comparison of the parameter-free LDA+DMFT(QMC) spectra of SrVO"3 (solid line) 
and CaVO"3 (dashed line) with experiments below and above the Fermi energy. Left panel: high- 
resolution PES for SrVO"3 (circles) and CaVO"3 (rectangles) Right panel: Is XAS for 
SrVO"3 (diamonds) and Cao.9Sro.1VO3 (triangles) Horizontal line: experimental subtraction 
of the background intensity. 
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